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ABSTRACT 

Early-type R stars and J stars are a special type of carbon star, having enhanced 
nitrogen ([N/Fc] « 0.5), lithium, a low ^^C/^^C ratio (< 15) and no s-element en- 
hancements. The merger of a helium white dwarf with a red giant is regarded to be 
a possible model for the origin of early-type R stars, but the details of nucleosyn- 
thesis are not clear. In this paper we investigate three possible channels for helium 
white-dwarf -I- red-giant mergers, and find that, amongst the three, only a high-mass 
helium white dwarf subducted into a low core-mass red giant can make an early-type 
R star. Nucleosynthesis of elements carbon, nitrogen, oxygen and lithium correspond 
well with the observations. Furthermore, we find that the J stars may represent a short 
and luminous stage in the evolution of an early-R star. 

Key words: stars: abundance, stars: white dwarfs, stars: evolution, stars: chemically 
peculiar, binaries: close 



1 INTRODUCTION 

Carbon stars are characterized by C/O > 1 in their atmo- 
spheres. They are generally classified as spectral types N, R 
and J. The N-type carbon stars are recognized as normal car- 
bon stars in which carbon is enhanced by third dredge-up 
during the asymptotic giant-branch (AGB) phase. R-type 
carbo n stars are further classified into two groups (|Shanel 
1 19281 ): early-R (hot) and late-R (cool). Late-R stars are 
similar to normal carbon stars (N-type), and the s-element 
over abundance is p ossibly produced during the AGB phase 
(IZamora et al ] |2009l '). Early-R stars are different from late-R 
stars; they have near solar metallicity and show enhanced 
nitrogen ([N/Fe] ~ 0.5), a low ^^C/"C ratio (< 15) and n o 
s-element enhancement (|Dominvlll98i : IZamora et al]|2009l ) . 
Also, the early-R stars are lithium-rich and none are found 
in binary systems. The J-type carbon stars are similar to 
early-R stars, but show a high luminosity, a smaller ratio 
ofJ£C/^^C and more lithium on the surface (|Abia fc IsernI 
I2OO0I ). The surface composition of N-type and late-R carbon 
stars ca n be reproduced by low-mass AGB nucleosynthesis 
models (|Zamora et al.|[2009l ). But the origins of the early-R 
stars and J stars are not clear, especially the nucleosynthesis 
responsible for their surface chemistries. The fact that, ex- 
cept for the wide binary BD+02°4338, all e arly-R stars are 
single stars, encouraged llzzard et al.l (|2007l 'l to argue for an 
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origin based on binary-star mergers, which was supported 
by binary-star population-synthesis calculations 

According to llzzard et al.l (|2007l ). a possible evolution 
channel for early-R stars is for a helium white dwarf (He 
WD) to merge with a red giant-branch (RGB) star. As the 
helium white dwarf comes into contact with the expand- 
ing red giant, a common envelope forms. Assuming spiral-in 
occurs faster than the envelope can be ejected, the white 
dwarf will merge with the helium core of the giant. The fi- 
nal outcome will be a single star having a core comprising 
the sum of the helium white dwarf and the re d-giant (RG) 
helium core surrounded by an H-rich mantle, llzzard et al.l 
(|2007h discuss the composition of the envelope, which seems 
to be able to reproduce a n enriched carbon surface. How- 
ever, [Piersa^r^r^al] (|2010l ) point out that the enhancement 
of carbon depends on whether efficient helium burning oc- 
curs during the merger process. To ana lyze the details of the 
merger process, iPiersanti et al.l (|2010l ) performed a three- 
dimensional smoothed-particle hydrodynamics (3D SPH) 
simulation and discussed the subsequent evolution. Accord- 
ing to their calculation, no efficient helium burning occurred 
which would dredge up carbon-enriched material to the sur- 
face. However, some high-mass helium-white-dwarf merger 
channels, i.e. where the whi te dwarf mass i s grea ter than 
0.2 M0, were not studied bv IPiersanti et all (|2010l '). 

In order to survey a wider range of models, we have 
simplified the SPH simulation into a ID stellar evolution 
calculation. A detailed nucleosynthesis analysis is included. 
Our result shows that the merger of a high-mass helium 
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white dwarf with a red giant could be a possible evohition 
channel for early-R stars and J stars. 



2 METHODS 

To inv estigate evolution following the merger of two white 
dwarfs, IZhang fc Jefenl (|2012bl ) considered three sets of as- 
sumptions concerning the rate of accretion onto the more 
massive star, representing different distributions of fast 
and slow accretion. To simulate the merger of a helium 
white dwarf with a red giant, f ollowing the results o f SPH 
simulations jYbon et ahl l2007l : iLoren- Aguilar et all |2009| : 
IPiersanti et aL 201Ct ). we consider only the case of fast ac- 
cretion (or a fast merger). 

Numerical simulations of stellar evolution are car- 
ried out using the stellar-evolution code Modules for 
Experiments in Stel lar Astrophysics (MESA version 4028) 
jPaxton et al.l l201ll ). For our experiments, it is difficult 
to control the required final mass of the helium white 
dwarf from a full binary star evolution calculation. Thus, 
an artificial method is adopted. We start with a zero-age 
main-sequence star of mass 2.0 Mq (metallicity Z = 0.02) 
and evolve it until the helium core reaches the required 
mass. Then the hydrogen envelope is completely removed 
to produce a naked helium core, essentially a pre-WD 
model. This is achieved artificially in MESA by switching 
off the nucleosynthesis and applying a high mass-loss rate 
until the required model is obtained. These models cannot 
ignite central helium burning and evolve straight to the 
white-dwarf (WD) track, ending up with a luminosity 
logL/L© = -2 (see I Zhang fc Jeffervl (|2012bl ') for details). 
Similarly, a model for the RG core is constructed in 
exactly the same way, except that we consider separately 
the cases where the RG core is either warm or cold. 
Here, a "warm" core represent a true model RG core 
which is only just degenerate (degeneracy parameter 
ri « 3(0.2 M0) - 10(0.3 M0)), whilst a cold core is much 
more degenerate (77 = 12(0.2 Mq) - 80(0.4 M©)). Our 
models do not include rotation; all mixing is complete and 
driven by convection (no semi-convection, no thermohaline 
mixing, no diffusion). We separate the subsequent simula- 
tion into three steps. 
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Step one: merge the helium white dwarf with t he he - 
lium core of the re d gian t. Foll owing lYoon et al.l (l20Q7h : 
iLoren-Aguilar et all (|2009l ') and IPiersanti et all (|201of K a 
fast accretion rate of 10* M0 yr~^ is adopted, implying that 
the merger takes a few minutes. 



Step two: Switch off time-dependent terms, including 
nucleosynthesis. Numerically, build up a hydrogen envelope 
on top of the merged helium core from step one. This 
is achieved by reversing the process used to remove the 
envelope, in this case by applying a negative mass-loss rate. 
The composition of the hydrogen envelope is defined to be 
identical to that of the RG progenitor. The structure of the 
post-merger model is thus a cool degenerate central helium 
core with a helium-rich hot shell surrounded by a hydrogen 
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Figure 1. The circles with crosses shows the three types of model 
computed. Dashed lines indic ate three zon e s defin ing possible evo- 
lution channels described bv llzzard et al. 



envelope. 



Step three: Switch on time-dependent terms and evolve 
the post-merger model. Stop evolution before the thermally- 
pulsing AGB phase starts. No stellar winds were included in 
any of the calculations. The post-merger models are defined 
to have a final mass of 2.0 M0; since the initial red giant 
has a mass of 2.0 Mq to which is added the WD mass, we 
assume that a fraction of the total has been ejected during 
the common-envelope phase, but have chosen this fraction 
arbitrarily for the time being. One calculation with a final 
mass of 1.5 Mq was made for comparison (§4.3.2). 



3 MODELS 



According to llzzard et al] (|2007l ). there are two main WD- 
RG merger channels: low-mass He WD -|- RG (R3a) and 
high- mass He WD -f RG (R3b). Channel R3a can be fur- 
ther divided into two groups: high core-mass RG and low 
core-mass RG. Thus, we set up three models for these three 
channels labelled Models A, B and C respectively (see Fig. [T] 

for details). 

In the IPiersanti et~all |20lO) SPH simulations, all the 
models are low-mass He WD + low core-mass RG mergers, 
and all the accretion is by mass transfer from the helium 
white dwarf to the RG core. However, one concern is that the 
He WD has a higher density than the RG core. The massive 
He WD may therefore subduct (or sink) into the low-mass 
helium core. To cover these possibilities, we have considered 
two separate merger types; i) accretion from the He WD 
onto the helium core and ii) accretion from the helium core 
onto the He WD, i.e. the He WD subducts into the center 
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Figure 2. Cartoon illustrating two possible outcomes of a wliitc- 
dwarf plus red-giant merger. 



of helium core; see Tableland Fig. [2] for details. As noted 
above, the composition of the hydrogen envelope is taken 
from a 2 Mq RG model which has a 0.2 Mq helium core, 
and the final masses are 2 M© . 



4 RESULTS AND DISCUSSION 

In both "accretion" and "subduction" models, non- 
degenerate helium is effectively added to the surface of a 
degenerate helium core at a high rate. The combination of 
gravitational contraction and compression by the outer lay- 
ers means that a hot shell develops at the WD/RG-core 
boundary, where the maximum temperature exceeds lO^K 
(Fig- El . Hydrogen burning starts from the helium/hydrogen 
boundary, where the released energy forces the star to ex- 
pand to large radius (Fig. [4]). The location of helium-ignition 
depends on the temperature and mass of the RG core and 
whether the WD is accreted or subducted. The merged 
stars are more luminous (or hotter) than normal red giants 
because of the larger core masses and, in some cases, an 
increased mean molecular weight in the envelope, due to 
dredge up of helium and carbon. 

In the cases that both the WD and the RG core are 
cold, a helium-burning flame starts burning from the merger 
boundary inward to the center of the core (Fig. (S)). As the 
helium-burning zone moves inwards, there follow a series 
of helium flashes, each subsequent flash decreasing in in- 
tensity, until the helium-burning shell (flame) reaches the 
center of the star. During the He-flash stages, luminosity 
first increases and then decreases. After the flame reaches 
the center of star, the star goes through a horizontal-branch 
or red-clump phase (core-helium plus shell-hydrogen burn- 
ing). After core- helium exhaustion, a helium-burning shell 
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Figure 3. The temperature profile around temperature maxi- 
mum for five models (labelled, see text and Table [TJ just after 
merger. In these models both the RG core and the white dwarf 
are cold. 
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Figure 4. The evolutionary tracks of all models. The different 
symbols with solid lines indicate the major evolutionary stage 
for each model (evolution is always toward higher luminosity). 
Symbols are separated by intervals of lO'^yrs. The first point on 
each track corresponds to lO^yrs after merger (marked by short 
horizontal ticks), the last is close to the start of the thermally- 
pulsing AGB; both can be identified in Figs.|6](A),[7](Bl),|8](B2), 
111! (CI), and 1121 (C2). The dashed line shows the evolutionary 
track of a normal 2 Mq star. 
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Table 1. Summary of the numerical experiments described in the text, showing the mass of the helium white dwarf, the core mass of 
the initial red giant, the assumed final mass of the merger product, and the type of white-dwarf-core interaction. 



Models 


WD mass (Mq) 


Core mass ( Mq) 


Final mass (Mq) 


Merger type 


Model A 


0.2 


0.2 


2.0 


accretion / subduction 


Model Bl 


0.2 


0.3 


2.0 


subduction 


Model B2 


0.2 


0.3 


2.0 


accretion 


Model CI 


0.4 


0.2 


2.0 


accretion 


Model C2 


0.4 


0.2 


2.0 


subduction 
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Figure 5. As Fig. [S] but for models where the RG core is warm 
and accretes a cold white dwarf. 



will again form and evolution will be similar to the AGB 
double-shell burning phase. 

Whilst our first calculations assumed a cold-degenerate 
model for the RG core, a more likely scenario is that the RG 
core should be warm at the time of merger (the core is heated 
by the hydrogen-burning shell). In practice, the pre-merger 
temperature of the accreted material is unimportant; fast 
accretion heats it to > 10*K very quickly. Hence for models 
where the WD is subducted into the RG core (Bl and C2 
below), the approximation is valid. In the remaining cases, 
we also computed models for accretion onto a warm RG core 
(Fig. [S]), and report the result for each case below. 

The majority of the test calculations do not show any 
change in the surface composition of the post-merger; where 
significant changes occur, we discuss the results in more de- 
tail. 



4.1 Model A: low-mass He WD + low-mass 
helium core 

Immediately after merger, hydrogen burning drives a fully 
convective envelope, whilst helium burning is ignited at 
0.2 Mq (the RG-core boundary) and moves inwards (Fig. [6]). 
^■^C is produced by helium burning through the 3a reac- 



Figure 6. The convection structure of Model A, a 0.2 Mq He WD 
accreted onto a 0.2 Mq RG core. The grey area shows the con- 
vection zones, the dashed line shows the hydrogen burning front 
and the solid lino shows the helium burning front. The bottom 
panel shows an enlargement around the phase of inward helium 
burning. 



tion during the hot accretion phase. During the inward 
helium-burning phase, there is no flash-driven convection 
zone which can reach the hydrogen-rich envelope. Thus, no 
burning ash (^^C) can be dredged up to the surface. In Fig.|6l 
the grey zone shows the convection zones, the dashed line 
shows the hydrogen-burning front and the solid line shows 
the helium burning front. Evolution takes 1.66 x 10* yr from 
merger to the start of the thermally-pulsing double-shell 
burning AGB. 

In the case of a WD accreted onto a warm RG core, 
helium ignition is nearly central, and more closely resem- 
bles the classical core-helium flash. There is no mixing of 
processed material to the sur face. 

These models confirm the lPiersanti et al] (|2010h conclu- 
sion that the merger of a low-mass He WD with a low-mass 
helium RG core cannot provide sufficient helium burning to 
bring ^^C to the envelope. 
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Figure 7. As Fig. [6l but siiowing tiic convection structure of 
Model Bl, a 0.2 M© He WD subducted under a 0.3 M© RG core. 
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Figure 8. As Fig. |6l but siiowing the convection structure of 
Model B2, a 0.2 M© He WD accreted onto a 0.3 M© RG core. 



4.2 Model B: low-mass He WD + high-mass 
helium core 

4.2.1 Model Bl (subduction) 

Model Bl (Fig.O is similar to model A, and no flash-driven 
convection zone reaches the hydrogen envelope. Thus, the 
evolution is similar except for helium-core mass (0.5 Mq 
compared with 0.4 Mq for case A) and evolution time (8.4 x 
10' yr). 



4.2.2 Model B2 (accretion) 

Because 0.3 M© is effectively accreted, the helium shell 
in model B2 is much hotter than in model Bl (where 
only 0.2 M0 is accreted), especially at the He/H boundary 
(Fig. As a consequence, hydrogen-burning at the bound- 
ary is very strong and drives a convection zone which mixes 
hydrogen downwards into the helium shell, and also moves 
the hydrogen- burning shell downwards by ~ 0.1 Mg, im- 
mediately after merger (t=0 yr) (Fig. [S]). This downwards 
mixing also dredges up CNO-cycled material which is mixed 
throughout the fully convective envelope, which thus be- 
comes C-poor and N-rich. It takes 7.6 x lO' yr from merger 
to the thermally-pulsing double-shell burning AGB. 

There are some globular cluster stars which are C- 
poor and N-rich which have not yet been fully explained 
jBrilev et al.ll2002l. |2004| : ICarretta et~ai]|2005l : ICohen et all 
I2OO5I : Bekki et al.ll200'if K Our model may account for some of 
these stars. Model B2 produces stars with surface [C/Fe] = 
-1.15, [N/Fe]= 0.88 and [Na/Fe]= 1.42, which compares 
with observed abundances of [C/Fe], [ N/Fe] and [Na/Fe] 



in the range of —0.5 to 0.5, —1 to 2 (ICohen et al 



Bekki et al 



Bekki et al 



2007 ) (Fig. [91) and 0.0 to 0.8 (|Yong et al 



20051: 



20031: 



20071 ). Model B2 is also helium rich, with Y 



0.4; some fraction of globular cluster stars have Y > 0.3 
(|Bekki et al.ll2007l ). 

The oxygen abundance is in the range observed in glob- 
ular cluster s (e.g. [0/Fe1=— . 10 in our model a nd —0.5 to 
0.5 in M15 (|Cohen et al.ll2005l : lBekki et al.ll2007^ ). However, 
our models are more C-poor than most C-poor globular- 
cluster stars, and also, our model produces lower-gravity 
and lower-temperature stars than observed in M15 (Fig. [TO}. 
However, we note that most of these cluster stars have lower 
metallicity than our models; it will be interesting to compute 
a low metallicity model in the future. 

In the case of a WD accreted onto a warm RG core, 
helium ignition is in the same place (the WD / RG-core 
boundary) but, because the core is non-degenerate and the 
maximum temperature at the WD / RG-core boundary is 
lower than that in the cold-core model (Fig. [5}, the shell 
flash is not strong enough to drive convection through the 
hydrogen shell, and there is no mixing of CNO-processed 
material to the surface. 



4.3 Model C: high-mass He WD + low-mass 
helium core 

4.3.1 Model CI (accretion) 

Model CI (Fig. [TTJ is similar to Models A and Bl and, 
again, no flash-driven convection reaches the hydrogen en- 
velope. Thus, evolution is similar except for a different core 
mass and evolution time (4.8 x lO' yr). The result for a WD 
accreted onto a warm RG core is also similar to that for case 
A, with no surface mixing. 
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Figure 9. [C/Fe] vs. [N/Fe] ratios for M odel B2 conipared 
with the globular stars of M15 (filled circles) llCohen et al.ll2005l ; 
iBekki et a,l1l2007lV The open circle with cross shows the model 
abundance of carbon and nitrogen during the major evolutionary 
stage. 
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Figure 10. Gravity - temperature diagram showing the r elatio n 
adopted for the g lobular stars of M15 by ICohen et al 
iBekki et al.l | |2007|) (thick solid line) and the evolutionary tracks 
for model B2 after merger (dashed line). The triangles denote 
intervals of lO'^ yrs during the major evolutionary stage (slow 
core- helium burning); other key times (since merger) are also in- 
dicated. 
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Figure 11. As Fig. |6l but showing the convection structure of 
Model CI, a 0.4 M© He WD accreted by a 0.2 M© RG core. 
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Figure 12. As Fig. |6l but showing the convection structure of 
Model C2, a 0.4 M© He WD subducted beneath a 0.2 Mq RG 
core. 



^.3.2. Model C2 (subduction) 

In model C2 (Fig. I12p . the helium shell is very hot, and en- 
hanced hydrogen-burning mixes ~ 0.2 Mq hydrogen into the 
helium shell. At the beginning of evolution, hydrogen burn- 
ing forms a fully convective envelope, and forces the star 
to expand. Immediately foUovying merger (and for approx- 
imately 150 years thereafter), the hydrogen-burning shell 
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converts available C to C and N, thus depleting C in 
the convective envelope and decreasing the ^^C /^'^C ratio 
by a factor of « 0.2 (the ^'^C /^'^C ratio in the envelope of the 
RGB progenitor is 26). Since the initial mixing is short-lived 
(50 years or four timesteps in our model), this factor is likely 
to be subject to several sources of uncertainty including the 
convection physics, numerics and the model parameters {e.g. 
masses and accretion rates). 

After this, the hydrogen-burning shell detaches from the 
convection zone. Meanwhile, the core is relatively cold (20 
MK), so the helium- burning shell is weak and relatively inac- 
tive {Lue/Lu ~ 10"*^), whilst the hydrogen shell burns out- 
wards through ~ 0.1 M0. After lO^yr, the He-shell temper- 
ature rises sufficiently for a helium-shell flash to occur, driv- 
ing a convection zone which mixes through to the hydrogen- 
rich envelope. Helium-burning continues in a series of flashes 
from the merger boundary to the center of the core. 

During the first two helium-shell fiashes, flash-driven 
convection connects the helium-burning shell to the convec- 
tive envelope and dredges freshly-produced ^^C ash from 
the helium-burning shell to the surface (the duration of 
these mixing events is short (< 200 yr) so that insufficient 
^^C(p,7)^^N(,/3+^)"C(p,7)^''N reactions occur as the car- 
bon transits the H-shell to affect the surface ^''C or ^^N 
abundances directly). However, the events are long enough 
compared with the convective turnover time, assuming that 
this is comparable with a normal red-giant envelope, that 
the envelope will be fully mixed. The evolution of the struc- 
ture is shown in Fig. 1131 flower panel), where the position of 
the H-burning shell is indicated by the maximum in nuclear 
energy production, which is very susceptible to local mix- 
ing events. Fig. [13] also shows the surface-abundance evolu- 
tion of several species over the same interval (upper panel), 
demonstrating the dredge- up of ^^C from 3q burning, en- 
hanced "N from CNO-cycle helium burning ([N/Fe]=0.69), 
■^^Ne from a captures on ^*N, and the depletion of ^^O, also 
associated with the CNO-cycle. The detailed evolution of 
the structures is likely to be sensitive to resolution and time 
step, as is the stability of the calculation; numerical prob- 
lems with forced changes to the time step prevented further 
exploration. 

The model takes 6.6 x 10^ yr to evolve from merger 
to the thermally-pulsing double-shell burning AGB; most of 
the time is spent as a red-clump star on the giant branch 
(squares in Fig. [TD) 

Fig. | 14 | compare s model C2 with 15 early-R stars from 
IZamora et al] (|2009l ') on a g — TbB diagram. The observed 
stars are satisfactorily close to the evolutionary tracks. How- 
ever, since most of our models (and models of other giants) 
are almost coincident on such a diagram, this is more of 
a consistency check than a strong constraint. Fig. [TS] com- 
pares the C/0 ratios and ^^C/" C ratios for Model C 2 with 
the early-R star measurements (|Zamora et al.l [20091 ). Dur- 
ing the first million years after merger, the model evolves 
through the region occupied by early-R stars. The precise 
trajectory is subject to the uncertainty in the initial deple- 
tion of ^^C discussed above, but the start and end points 
should be similar; it is conjectured that the allowable tra- 
jectory envelope should bracket the loci of observed early-R 
stars. The models generate a similar abundance of C and O 
to that observed. Furthermore, we find that ^^Ne and ^^Na 
are enriched to give [^^Ne/Fe] = 1.98 (^°Ne/^^Ne= 0.17) and 
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Figure 13. Bottom: as Fig. 1121 but expanded around the the 
first two helium-shell flashes of Model C2. Top: the evolution of 
the surface abundance of ^^C, ^*N, ^''O and ^^Ne over the same 
interval. 



-1 



- 



in 

E 
o 
3 

O 2 



A'' 

■ — 1 Myr 




_J_ 



3.8 



3.7 



3.6 

log(Tea/K) 



3.5 



3.4 



Figure 14. Gravity - temperature diagr am showing the early- 
type R stars (circles) IZamora et al.l | |2009| ) and the evolutionary 
tracks for model C2 (dashed line). Squares denote intervals of 
10^ yrs during the major evolutionary stage. The cross shows the 
mean errors for the R-star measurements. 
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Figure 15. The post-merger evolution of the C/O and ^^C/^^C 
ratios for Model C2 compared with the early-R stars of 
IZamora et al] | |2009|') . The star HIP 58 786 has I2c/l3c=70 and 
is not shown on this diagram. The dashed line shows the model 
evolution of C/O and ^^C/^^C ratios after merger. The circle 
with cross shows the ratios of the C/O and ^^C/^^C during most 
of this evolution. The cross shows the mean errors for the R-star 
measurements. 



Figure 16. Lithium versu s ^^C /^^C. Circles show the early- 
R stars from [Zamora et al. ] l l200a) . Dots show the J stars from 
lAbia fc Isernl l l200d) ." The dashed line shows the Li/H and 
12(~<^13q ratios during the evolution of model C2 immediately af- 
ter merger. The circle with cross shows the Li/H and ^^C/^^C ra- 
tios during the main core-helium burning phase. The cross shows 
the mean errors for the R-star measurements. 



[Na/Fe] = 1.46. It is not known whether early-R stars are en- 
riched in these species. Although the ^■^Ne(a,n)^^Mg is a 
putati ve neutron sour ce for further (s-process) nucleosyn- 
thesis (|Cameronl 1 196Cf ) . conditions in our models (Tmax ~ 
1.4 X lO^K with p « 3.8 X lO^g/cm^ in the helium sheU) do 
not satisfy the neccessary criteria (T > 2.2 — 3 .5 x lO'^K an d 
p > 1 - 3 X lO^'g/cm^) jRaiteri et al.lll99ll : iMevej Il994h . 
Moreover, s-process nucleosynthesis is not included in the 
version of MESA used in these calculations. 

For Model C2 with a final mass of 1.5 M© the evolution- 
ary tracks are similar to those with a final mass of 2.0 M©, 
except that the final C/O ratio increases to 6.1 (instead of 
4.2). This is because both models have the same core mass, 
and dredge up the same mass of carbon, which is conse- 
quently more concentrated in the lower-mass envelope. 



4.4 J stars and Lithium 

Early-R and J stars have a notably large overabundance of 
lithium, which has so far been difficult to ex plain. To explain 
the ov erabundance of lithium in RCrB stars. iLongland et al.l 
(|2012l ) showed that lithium can be produced by the merging 
of a helium white dwarf with a carbon-oxygen white dwarf if 
their che mical composition i s rich i n ^He from the previous 
evolution. IZhang fc Jefenl (|2012al ) confirmed this in their 
He+He WD merger simulation. This model requires enough 
■^He to be left in the white dwarf after the end of main- 
sequence evolution and a hot enough zone to form during 
the merger. In our simulation, we obtained a mass fraction 



of 5.6 X 10~* ■'He in the RG envelope before merger. We 
introduced this mass fraction into the envelope of the post- 
merger models. The hot He-burning shell drives convection 
which mixes ^He down into the He-burning shell where it 
burns with ''He to produce ^Li pHe('*He, 7)''Be(e" , i/)''Li]. 

For the new born ^Li to survive, the ^Be must be 
transferred to a cooler region rapidly before the decay 
occurs. Thus, in a manner similar to that described by 
ICameron fc Fowled (Il97iri . the Li may or may not ap- 
pear on the surface depending on whether sufficient con- 
vection occurs (see Fig. I12|l . In our models. Models B2 and 
C2 both produce lithium. However, in Model B2 the new- 
born lithium is almost totally destroyed immediately by the 
reaction ^Li(p, ''He)''He. In Model C2, the star retains its 
high lithium abundance because the mixing is very effective. 
While the star goes through its high- luminosity phase, its 
surface reaches a peak lithium abundance while the ^^C/^'^C 
ratio is still low (Fig. I16p . Thus, J stars may represent just 
a special stage dur ing the formation of early-R stars, as al- 
ready hinted at bv lAbia fc Isernl (|2000t ). 

This stage, where ^^C/^^C is very low, takes around 1 
Myr, compared with 60 Myrs for the lifetime of Model C2. 
If C2 represents early-R stars, then J stars should constitute 
PS 1/60 of the total (Fig. fTT)) . 

Although the statistical evidence is not very strong, 
there are beli eved to be ar ound 10 times more early-R stars 
than N stars (|Blancdll965l ). If the J stars constitute « 1/60 
of early-R stars, then the J:N ratio would be ~ 15% w hich 
is the observed value jAbia fc Isernll2000l : lBlancdll965l 'l. 
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WD subducting into a low-core-mass RG star (Model C2) 
can produce early-R stars. However, nucleosynthesis of ele- 
ments C, N, O, and Li for this model match the observations 
well. Moreover, J-stars may represent a short stage in this 
evolution. The fact th at only Mode l C2 c an produce early-R 
stars may explain whv llzzard et al] (|2007l ) produce ten times 
as many early-R stars as required. Our study shows that two 
of their proposed merger channels do not produce the correct 
surface composition. A more detailed population-synthesis 
study should be carried out to distinguish the boundaries of 
each of the evolution channels. 

Our study has also shown that a low-mass He WD ac- 
creted onto a high-mass RG core may account for some C- 
poor and N-rich globular cluster stars. Lithium enrichment 
in red giants has previously been explained by one or more 
of rotationa l mixing, mass loss, pl anet accretion and deep 
circulation (jCharbonnel fc Balachan dran 2000). In this pa- 
per we have shown that a He WD + RG merger may be an 
additional and important channel for producing lithium-rich 
red giants. 



Figure 17. The evolutionary track of Model C2. Squares repre- 
sent intervals of lO^yrs. The double-ended arrows show the lumi- 
nosity ranges for model C2 over the time intervals indicated on 
the figure. These correspond approximately with the observed lu- 
minosity ranges of early-R and J stars. The dot shows HD 100764. 
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This argument assumes that J-stars, like early-R stars, 
are all single. This may not be true. 

Approximatel y 10-15% of J stars are believed to have a 
silicate dust shell (|Llovd Evanslll99ll ). Some consider such 
shells to be circumbinary disks. Although there is evidence 
for disk asymmetry, no J -star companion, o r J-star orbit has 
so far been established (|Deroo et al.ll2007t l. If we consider 
the J-star phase of Model C2 to represent the immediate 
post-merger phase, it is possible that the observe silicate 
dust shell may form from the ejected part of the envelope. 
Since energy considerations demand that some fraction of 
the envelope should be ejected during merger, an asymmet- 
ric shell (or disk) would be an obvious geometry for the 
ejecta. If a silicate dust shell forms during merger and has 
a lifetime « 10^ years, then we might expect a fraction of 
10-15% to show such shells, since the J stars are younger 
than the majority of early-R stars and have a lifetime of 
~ IMyr. We still need to explain why HD 10 0764, an early- 
R star fFig. | 17J). also shows a dust shell jParthasarathvl 
I1991I : [sioan et al.ll2007l V However it has quite a low ratio of 
i2Qyi3(-i_^^ and show s emissions from polycyclic aromatic 
hydrocarbons (PAHs) |sioanelai][2003). We therefore sug- 
gest that HD100764 might be a very recent merger currently 
evolving towards the J-star phase. 



5 CONCLUSION 

Mergers of helium white dwarfs with red giants have been 
considered as a possible model for the origin of early-R stars. 
In this paper we have investigated three possible channels 
for He WD + RG mergers, and find that only a high-mass He 
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